Phase behavior of two-component lipid membranes: theory and experiments 
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The structure of the ripple phase of phospholipid membranes remains poorly understood in spite 
of a large number of theoretical studies, with many experimentally established structural features of 
this phase unaccounted for. In this article we present a phenomenological theory of phase transitions 
in single- and two-component achiral lipid membranes in terms of two coupled order parameters - 
a scalar order parameter describing lipid chain melting, and a vector order parameter describing 
the tilt of the hydrocarbon chains below the chain-melting transition. This model reproduces all the 
salient structural features of the ripple phase, providing a unified description of the phase diagram 
and microstructure. In addition, it predicts a variant of this phase which does not seem to have 
been experimentally observed. Using this model we have calculated generic phase diagrams of 
two-component membranes. We have also determined the phase diagram of a two-component lipid 
membrane from x-ray diffraction studies on aligned multilayers. This phase diagram is found to be 
in good agreement with that calculated from the model. 

PACS numbers: 87.16.D-,61.30.Dk 
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I. INTRODUCTION 

Phospholipids are amphiphilic molecules with a hy- 
drophilic head group and one or more hydrophobic hydro- 
carbon chains, and are the major constituents of biomem- 
branes. They self-assemble in water to form a variety of 
thermodynamic phases, depending upon the temperature 
and water content [l| . The most common of which is the 
lamellar phase, a periodic structure consisting of lipid 
bilayers separated by layers of water. 

Above the chain melting transition (also known as the 
main-transition) temperature T m , phospholipids form 
the fluid L a phase, where the hydrocarbon chains are dis- 
ordered with many trans-gauche isomerizations in their 
C-C bonds. Positional order in the plane of the bilayer is 
liquid-like in this phase. The lower temperature gel (Lp 
or Lpi) phase is characterized by flat bilayers with the 
hydrocarbon chains predominantly in the fully-stretched 
all-trans conformation. In the Lpi phase the chains are 
tilted relative to the bilayer normal, whereas in the Lp 
phase there is no such tilt. The presence of a tilt depends 
on the relative cross-sectional areas of the head group 
and chains of the lipid species @, Q • For example, phos- 
phatidylcholines display a tilt angle of about 30° while 
phosphatidylethanolamincs with relatively smaller head 
groups do not show any tilt. 

Phospholipids such as phosphatidylcholines (PC) ex- 
hibit a modulated phase between the L a and the Lpt 
phases at high water content [l|, Ij-Q ■ This is known as 
the Ppi or the ripple phase, and is characterized by a 
one-dimensional periodic height modulation of the bilay- 
ers. In this phase the hydrocarbon chains are partially 
ordered and have a non-zero tilt relative to the bilayer 
normal. The first order Lpi - L a transition is well un- 
derstood in terms of the melting of some of the degrees 
of freedom within the hydrocarbon chains [7]. However, 
the first-order Lpi - Pp* transition (often called the pre- 



transition), which occurs at a temperature slightly below 
the main transition, is much less understood. 

The structural features of the ripple phase have been 
extensively studied using a variet y o f techniques such 
as, x-ray diffraction (XRD) [l], l8-flll]. neutron diffrac- 
tion [lj| , freeze- fracture electron microscopy [H Il3ril6| , 
and scanning-tunneling microscopy (fj, \\J\ . These exper- 
iments have revealed the existence of both symmetric as 
well as asymmetric ripples. The asymmetric structure 
is, however, thermodynamically more stable and is com- 
monly reported in almost all x-ray studies on PC sys- 
tems. The wavelen gth of the ripples typically lies in the 
range of 12-16 nm [18| and the peak-to-peak amplitude 
is about 2 nm. X-ray diffraction patterns from this phase 
can be indexed on a two dimensional oblique lattice and 
the rippled bilayers lack a mirror plane normal to the 
rippling direction. Calculations of the electron density 
distributions [ll|, 120-23] confirm the basic asymmetry of 
the ripples, and indicate that the bilayer thickness is ap- 
preciably different in the two arms of the ripple. The ab- 
sence of a mirror plane normal to the ripple wave vector 
in the asymmetric ripple phase can be attributed to the 
presence of asymmetry in either the shape (i.e., one arm 
being longer than the other) or bilayer thickness (with 
a higher bilayer thickness in one of the arms) or both. 
However, in almost all the systems studied so far both 
these features are present simultaneously. 

Several theoretical models have been proposed to de- 
scribe the ripple phase [2443 lj , but none of them account 
for all the experimental observations, namely, (i) occur- 
rence of the Ppi phase between L a and Lp> phases, sep- 
arated by two first-order transitions; (ii) unequal bilayer 
thickness in the two arms of the ripple; and (iii) unequal 
lengths of the two arms. 

Macroscopic continuum theories, which treat the bi- 
layer as a continuous membrane, can broadly be classi- 
fied into two main categories depending on the nature of 



the order parameter that is chosen to describe the ripple 
phase. One set of theories considers the order parameter 
to be a scalar, such as the membrane thickness as in the 
models by Falkovitz et.al [25fl. Marder et.al [26| and by 
Goldstein and Leibler (GL) [28|, [29| or the configuration 
of the hydrocarbon chains (ratio of cis to trans bonds in 
the chains) |32|]. Due to the scalar nature of the order 
parameter, these models cannot distinguish the different 
symmetries of the L a , Lpi and Ppi phases. Further, in 
these models the ripple phase is characterized by a bi- 
layer thickness modulation and not a height modulation. 
The other set of continuum theories considers the order 
parameter to be a vector, such as the tilt of the hydro- 
carbon chains with respect to the local layer normal, as 
in the Lubensky-MacKintosh (LM) [30] and the Seifert- 
Shillcock-Nelson models [33j |. 

Although models involving vectorial order parameter 
can predict the various phases with different in-plane 
symmetries, they have their own limitations. For exam- 
ple, the LM-model cannot account for the occurrence of 
the ripple phase in bilayers composed of achiral molecules 
and speculated that the the origins of the experimen- 
tally observed asymmetric ripples lay in the chirality of 
lipid molecules. However, subsequent experiments using 
racemic (achiral) mixtures showed this was not the case 

MM- 

In this paper we present a phenomenological Landau 
model to describe the ripple phase, that unifies the two 
distinct approaches (those using scalar and vector order 
parameters) used in the earlier models. In addition to 
the coupling between the curvature of the bilayer to the 
divergence of the molecular tilt [30( , we bring out the im- 
portance of the coupling of variation in the bilayer thick- 
ness to molecular tilt and show that the (i) asymmetry 
of the ripple phase arises from (a) the difference in the 
bilayer thickness in the two arms of the ripple, (b) un- 
equal length of the two arms, (ii) existence of a nonzero 
mean tilt of the hydrocarbon chains of the lipid molecules 
along the direction of ripple. In section II we present 
the model for single-component achiral membranes. The 
main results of this section has already been reported in 
a short communication 36] . In section III we extend the 
model to two-component membranes and use it to calcu- 
late the phase diagrams for various binary lipid mixtures. 
In section IV we present the phase diagram of dipalmitoyl 
phosphatidylcholine (DPPC) - dimyristoyl phosphatidyl- 
choline (DMPC) mixtures determined from small-angle 
x-ray scattering studies of aligned samples. The experi- 
mental phase diagram is found to be in good agreement 
with the one calculated using our model. 



II. SINGLE-COMPONENT ACHIRAL 
MEMBRANES 

In this section, we develop a phenomenological model 
for describing the asymmetric ripple phase in achiral 
membranes. To simplify matters, we ignore interbilayer 



interactions and choose to describe the ripple phase in 
terms of (i) a scalar order parameter ip, which represents 
the relative bilayer thickness and (ii) a vector order pa- 
rameter m, which is the projection of the molecular axis 
in the plane tangent to the bilayer surface (fig.©) [3Cj . 
We consider three distinct contributions to the total free 
energy density: 

a. The stretch free energy density f st : 

The main transition is characterized by the freezing 
of the lateral motion as well as conformational order- 
ing of the hydrocarbon chains. Many experiments have 
revealed that at the main transition there is a discontin- 
uous jump in the membrane thickness 37[. Therefore, 



following GL, we describe the main transition as one in- 
volving a change in the bilayer thickness and introduce a 
dimensionless scalar order parameter ip, such that: 



v> 
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where 6(r) is the actual membrane thickness, So is the 
constant thickness of the membrane in the L a phase and 
r denotes the lateral position within the bilayer plane. 
No specific assumption is made about the microscopic 
state of the lipids in choosing ip as the order parameter; 
we can think of ip as a hybrid coarse grained variable 
that summarizes changes in various types of degrees of 
freedom, like for example, conformational changes of the 
hydrocarbon chains, ip is taken to be positive for T < 
T m due to the stretching of the chains. This is valid in 
general, even if the chains are tilted below T m . 

The stretching free energy per unit area, for an isolated 
lipid bilayer membrane, can therefore be written as |38| : 



fst =~ a 2 i? + ~ a 3 V 3 + \ «4 ^ + \C (V^) 2 



\d(v^) 2 



\e (V^) 4 



(2) 



The explicit temperature dependence is assumed to re- 
side solely in the coefficient of ip 2 : 



do 



= a' 2 (T- T*) 



T* being a reference temperature (which in the absence of 
the cubic term in f st is the critical temperature). Since, 
Lp> — » L a phase transition is known to be first order, the 
coefficient 03 is taken to be negative. In the presence of 
the cubic term, the critical temperature T* is preempted 
by a first order melting transition at, 



T — T* 
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The coefficient C can either be positive or negative but 
D , E and 0,4 are always positive to ensure the stability 
of this free energy expansion. With C > 0, the equi- 
librium phases are always spatially homogeneous; either 
as L a or Lp (L«')- However, with C < modulated 



phases are possible with some characteristic wave vector 
q Q . To motivate how C can be negative, it is convenient 
to introduce an auxiliary scalar variable p. Constructing 
effective cylinders enclosing the head group and the chain 
with cross-sectional areas ah and a t , respectively, we de- 
fine p = (ah — at)/ (ah + at). Consider a bilayer membrane 
whose thickness ip increases steadily over a length scale 
£. This variation in thickness produces a strain at the 
molecular level. This can be accommodated both by a 
variation in the tilt m, and a variation in the mismatch 
p over this scale. The latter coupling to lowest order 
is of the form Vtp ■ Vp. This term, consistent with the 
symmetries of the bilayer, will generically lead to C < 0. 

b. The tilt free energy density ftut '■ 




FIG. 1: The unit vector n represents the orientation of 
the long axis of the lipid molecules relative to N, the bi- 
layer normal. The projection of n on the bilayer plane 
m (= n — (N ■ n) N) is the order parameter. m x and m y 
are the components of m along the two orthogonal directions 
in the plane of the bilayer. 

The tilt free energy density can then be written in 
terms of m as, 

ftut =\bi |m| 2 + h> A |m| 4 + Ti (V ■ m) 2 

+ r 2 (v 2 m) 2 + r 3 (v • m) 4 + r^ H 2 
+ r 5 (m • w-) 2 + r 6 (m x w) 2 



stretching free energy density and the tilt free energy 
density, respectively. 

As 62 is positive in our model, the term involving the 
coefficient T4, which couples the molecular tilt to the 
change in the bilayer thickness, is the most crucial one 
for the appearance of the asymmetric ripple phase. When 
T4 > 0, the stable phase below T m is Lp with |m|=0. On 
the other hand, tilted phases can form if T± < 0. The 
succeeding two terms take into account the anisotropy of 
the tilted bilayer. The next two terms represent higher 
order couplings between modulations in ip and in m, al- 
lowed by the symmetry of the system. These higher order 
coupling terms are kept merely for consistency and are 
not crucial for obtaining the optimal phases. 

c. The curvature energy density f c : 

The curvature energy density of the bilayer can be 
written as 1301; 



fc 



1 



(V 2 /i) 2 - 7 (V 2 /i) (V-m) 



(5) 



where h(x, y) is the height of the bilayer relative to a flat 
reference plane, k is the bending rigidity of the mem- 
brane, and the coefficient 7 couples the mean curvature 
to splay in m. 

The total free energy density is given by the sum of the 
above three contributions : 



/ — fst + ftut + fc 



(6) 



The equilibrium height profile of the bilayer h(x,y) is 
related to the tilt m via the Euler-Lagrange Equation, 



V 2 h = 1 (V • m) 



(7) 



Eliminating h from the free energy density / leads to the 
effective energy density / e // with a reduced splay elastic 
constant. 



T l =T l 
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+ T 7 N\) (V • m) 2 + T 8 (VVO 2 (V • m) 2 (4) Thus ' the effective free energy density is given by, 



First and second terms in the above equation represent 
the coupling of the molecular tilt vector m to the bilayer 
thickness and modulations in the thickness, respectively. 
The third term corresponds to the splay deformation of 
the m field and represent the tilt elasticity in the tilted 
phase. Since the order parameter m is a vector, cubic 
terms in m are not allowed. 

As shown by Jacobs et. al [39| m the context of a one 
dimensional model with a scalar order parameter, a term 
in the fourth power of the gradient of the order param- 
eter is necessary to stabilize a modulated phase with a 
non zero mean value of the order parameter. Hence we 
include (Vi/>) and (V ■ m) terms in our expression for 



feff = \ a 2 ^ + I a 3 V 3 + \a A ^ + l - C (V^f 



+ ■ l -D(V 2 4,f+ -E(V^t+ -U 2 |mr 



+ -b 4 |m| 4 + f 1 (V • m) 2 + T 2 (V 2 m) 2 

+ r 3 (v • m) 4 + r 4 v> |m| 2 + T 5 (m • v</>) 2 

+ T 6 (m x Vi/0 2 + r 7 (vV) (V ■ m) 2 

+ r 8 (wo 2 (V • m) 2 (8) 

To determine the mean field phase diagram we assume 
that the modulated phase is characterized by the most 



dominant wave vector q, and choose the following ansatz 
for ip and m, 

ip = ipo + ipisin(qx) 
m x = rriQ X + m\ x cos (qx) + m 2x sin (qx) 

+ n%3 x cos (2qx) + rriAx sin (2qx) 
m y = m 0y (9) 

with the amplitudes {ipi,m,i} as variational parameters. 
We retain Fourier components of m x to second order to 
account for the ripple asymmetry; the experimentally ob- 
served asymmetric ripple profile is obtained if m 3x ^= 0. 
In principle, one could have done a numerical variational 
calculation with many more Fourier components. How- 
ever, our simplified ansatz is sufficient to recover all the 
qualitative features of the experimentally observed asym- 
metric ripple phase. 

2D modulations of the membrane [40] was not included 
in our ansatz (eqn.[9]), since earlier studies [22|, [3fJ, |31j 
based on a similar free energy density for the tilt and 
curvature alone, show that the presence of a mean tilt 
suppresses two dimensional height modulations. The ad- 
ditional ip dependent terms in the free energy density do 
not alter this conclusion. Hence we confine our atten- 
tion to one-dimensionally modulated ripples. Note that 
spatial modulations in m y are neglected as we do not 
keep terms proportional to (V x m) in eqn.Q. Varia- 
tion of rriy with x represents non-parallel orientations of 
the lipid chains, which are not favored by van der Waal's 
interactions between the chains. Such terms are impor- 
tant in chiral systems, where they can lead to structures 
with non-zero winding numbers [30| . In the absence of 
chirality, these terms can, however, be neglected. 

The average free energy is obtained by integrating the 
F e ff over one spatial period, i.e, 



(feff) 



Jo' feffdx 

2-n- 

Jo* dx 



(10) 



A. Results and discussion 

The phase diagram in the C-T plane is obtained from 
the numerical minimization of the effective free energy 
(feff)- Fig. ([2]) represents an example of such a phase 
diagram. It is calculated for a choice of parameter val- 
ues, which reproduce closely the main- and pre-transition 
temperatures of DPPC at C = - 4.84xl0" 19 J. These 
values (except for Ti to Tq) are taken from ref [29| and 
ref [2a]. It is seen from fig.© that there are three dis- 
tinct regions in the phase diagram corresponding to three 
different phases: 



L a phase : 



V>o = 0, Vi = 0, m 0x = 0, m 0y = 
m lx = 0, m 2x = 0, m 3x = 0, m ix = 
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FIG. 2: Phase diagram in the T — C plane calculated from 
the model. a' 2 = 159.42 k B , T* = 260.0 K. Values of the other 
coefficients in units of ksT* are: 0,3 — -306.5, 04 = 613.15, 
62 = 0.2, b 4 = 200.0, D = 557.41, E = 600.0, Ih = 0.010, T 2 

= 1.80, r 3 = 500.0, r 4 = -3.0, r 5 = -20.0, r 6 = -20.0, r 7 = 

-50.0, Tg = -25.0. Both the main-transition (L 



Ppi) and pre-transition (Pp 



-'P' 



are first order. 



Lpt; L 



La' phase : 



Ppi phase : 



ipo ^ 0, Vu = 0, m 0x ^ 0, moy / 



TOi a 



0, m 2x = 0, m 3x = 0, mix = 



bo ^ 0, Vi ¥" 0, m ax ^ 0, m 0y = 
m x = 0, m 2x ^ 0, m 3x ^ 0, m ix = 



The first-order transition lines which separate these 
three phases meet at a Lifshitz point at C^ v — —3.43 x 
10" 19 J and T = 38.75°C. When \C\ < \C Lp \, the first 
order transition line is parallel to the C axis and occurs at 
T = 38.75°C. As long as \C\ < \C Lp \ the Lp phase melts 
directly to the L a phase at T m . But for \C\ > |Cl p | the 
La' phase first melts into the Ppi phase, which in turn 
transforms into the L a phase at higher temperatures. It 
is also clear from the figure that, as |C| increases, the 
region of the Ppi phase expands at the expense of the 
Lpi and the L a phases. 

It would seem that our free energy functional contains 
a large number of parameters, which we need to fine tune 
to obtain the asymmetric ripple phase. We would like to 
emphasize that asymmetric ripples are a fairly generic 
and robust feature of our model whenever T^ < 0. The 
other r^'s only make the shape anisotropy in the height 
profile more pronounced when they are negative, and 
are not crucial for the existence of the asymmetric rip- 
ple phase. Further, x-ray diffraction experiments identify 
the asymmetric ripple phase on the basis of its oblique 
unit cell. Such a unit cell can arise either from a strongly 
asymmetric height profile or from a weakly asymmetric 
height profile as long as the thickness is different in the 
two arms of the ripple [22] . This feature is exactly what 
is captured within our theory (Figs. 4(a) and|4(b)). 
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FIG. 3: Dependence of the phase behavior on Ty. All lines in 
the phase diagram correspond to first-order transitions. C = 
5.74X10 -19 J, and all other parameters as in fig[2] 



In addition to C and 1^, the phase behavior depends 
crucially on IV Depending on the value of Ty three dis- 
tinct modulated phases are found, with differing spatial 
variations of the order parameters (ng[3J). For T 7 < 0, 
the P/3' phase is stable, where the height profile is asym- 
metric and resembles very closely those seen in experi- 
ments (fig 4(a)) [17h19J . In this phase m\ x = 0, mix. ^ 
0, m 3x ^ 0, m,4 X = 0. -01 is almost 7r/2 out of phase 
with h, so that it is positive (negative) along the longer 
(shorter) arm of the ripple, resulting in different bilayer 
thicknesses in the two arms, again in agreement with 
experimental observations [f8|, fL3] . W hen Ty is weakly 
positive, the Pg, is obtained (fig |4(b)| , where the higher 
Fourier components of m x are absent (m\ x — 0, m 2a; ^ 
0, ra^x = 0, m4 X = 0) and hence the height profile is 
symmetric. As in the Ppi phase the bilayer thickness 
is different in the two arms. This profile resembles the 
structure reported in ref. [22], based on the diffraction 
data of ref. [1J. For higher values of Ty the P 2 p> is found 
(fig|4(c)|>, where m lx ^ 0, m 2x = 0, m 3x = 0, m 4x ^ 0. 



The height profile is symmetric and is almost in phase 
with ipi. As a result, the bilayer thickness is modulated 
within each arm of the ripple. This novel structure does 
not seem to have been observed in any experiments till 
now. However, since the electron density map of the rip- 
ple phase has been calculated only in a few instances, 
the possible occurrence of this structure in some systems 
cannot be ruled out. 

The present model predicts the chain tilt to be high 
in the thicker arm and negligible in the thinner one. 
There are in principle at least two chain conformations 
consistent with this thinner arm of the bilayer - one is 
by having disordered chains, the other by chain inter- 
digitation. Since in our model, the high temperature 
phase was identified as L a , it is natural to populate the 
thinner arm with disordered chains. We therefore pro- 
pose that in the asymmetric ripple, the thicker arm is 
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FIG. 4: Spatial variation of if), m x and h in the different 
ripple phases obta ined for |(a)| Yy = —150 ksT* , [(b)] F 7 = 
100 k B T* and [(c)] T 7 = +150 k B T* . T = 310.50 K and all 



other parameters are the same as those used in figE] 



made up of tilted gel-like domains and the thinner one 
is fluid- like with negligible chain tilt. This conclusion is 
partially supported by spectroscopy and diffusion experi- 



ments, which indicate a significant fraction of disordered 
chains in the ripple phase [4l|,[42|. However, it must be 
acknowledged that there is no direct experimental data 
on the spatial distribution of the disordered chains within 
the rippled bilayer. Electron density maps calculated 
from diffraction data do not provide any direct informa- 
tion about chain ordering. In contrast, recent computer 
simulation studies of the ripple phase using a variety of 
model lipids at different levels of coarse-graining J43H4~5|. 
find an asymmetric ripple structure consistent with the 
height and thickness modulations seen in the electron 
density maps, with the thicker arm consisting of gel-like 
domains, and the thinner one made up of ordered and in- 
tcrdigitated chains. The cross-over region between these 
two contains a large fraction of disordered chains. In the 
gel-like domains chains are tilted with respect to the lo- 
cal layer normal, whereas chain tilt vanishes in the thin- 
ner arm due to interdigitation. Thus the height, thick- 
ness and tilt modulations found in these simulations are 
consistent with the predictions of the present model, al- 
though the thinner arm consists of interdigitated chains 
instead disordered chains as proposed here. Symmet- 
ric ripples consisting of alternating gel-like and fluid-like 
domains have also been found in some simulations [46j . 
These ripples are characterized by only a thickness mod- 
ulation and no height modulation. Stacks of such bilay- 
ers will result in a centered rectangular lattice and do 
not seem to have been observed in any pure lipid-water 
system. 



III. BINARY LIPID MIXTURES 



system there is a contribution f m ix to the total free en- 
ergy, corresponding to the free energy of mixing. f m ix is 
given by the sum of the entropy of mixing and enthalpy, 
which within the Bragg- Williams (mean-field) approxi- 
mation can be written as, 



f„ 



k B T[0log(0) + {l-0)log(l 



+ ij0(l-< 



(11) 
where, k B is the Boltzmann constant, and J > is an 
attractive interaction parameter that enhances lipid-lipid 
demixing. 

Since the important properties of the lipids are mainly 
reflected in the model parameters T* , C and r 4 , we as- 
sume their values for a mixture to be the average of that 
for the two pure lipids weighted by their concentration. 



C{<£) = 


<l>c A + (\-<p)c B 


(12) 


r 4 (0) = 


4> v iA + (l - 4>) r 4s 


(13) 


T*(4>) = 


<t>T* A + {l-4,)T B 


(14) 



where, the subscripts A and B refer to the lipids A and B, 
respectively, and </> is the mole fraction of lipid A in the 
mixture. All other parameters are assumed to be same 
for both the lipids. 

Combining eqn.(j8]) and eqn.lfTTj). we obtain the total free 
energy density f total as, 



Jtotal Jmix > Jeff 



(15) 



After minimizing Eqn.(|15[) with respect to both tp and m, 
the two-phase coexistence in the (T, </>) plane is obtained 
using Maxwell's construction. 



Biomcmbranes are multicomponent, including lipid 
species that show liquid-ordered behavior, and are likely 
to have thickness variations [47]]. The first step in un- 
derstanding such complex systems is to study binary 
mixtures of lipids which differ in their head groups and 
chains. Therefore, in this section, we extend the model 
presented in the previous section for single component 
membranes to membranes consisting of two lipid species. 

We specifically consider two different cases. In case 
1, only one of the lipids exhibits the ripple phase. Since 
the occurrence of the ripple phase is closely related to the 
formation of a tilted gel phase at lower temperatures, one 
of the lipids is assumed to show the phase sequence L a — > 
P/3> —¥ Lp> , and the other to show L a — >■ Lp on cooling. 
This case corresponds to two species differing only in 
their head groups. In case 2, both the lipids exhibit the 
ripple phase, and corresponds to the situation where the 
two lipids differ only in the length of their hydrocarbon 
chains. We employ an approach quite similar to that 
used in ref 38[ and ignore any lipid exchange with the 
surrounding solvent. Owing to the fact that the lipids 
constituting the mixture can have different chain lengths, 
degrees of saturation, or head groups, the main transition 
temperature is taken to be different for the two lipids. 

Due to the presence of an additional component in the 



A. Results and discussion 

1. Tilt-Non-Tilt Mixtures 

We now calculate the phase diagram of binary 
lipid mixtures, where only one of the constituents 
exhibits the rippled phase, such as mixtures of di- 
palmitoyl phosphatidylcholine (DPPC) and dipalmitoyl 
phosphatidylethanolamine (DPPE). The phosphatidyl- 
cholines, which have the —N(CH^)^ moiety in their po- 
lar head groups, typically exhibit a tilted Lp> phase. On 
the other hand, the phosphatidylethanolamines which 
have the smaller —(NH^)^ moiety instead, shows the 
untilted Lp phase and do not exhibit the ripple phase. 

From fig ([2j one finds that only when the value of \C\ 
is greater than |Cl p |, a spatially modulated phase occurs 
between the high and low temperature phases. The sign 
of T4 is also important in determining the existence of a 
tilted phase at lower temperatures. If T^ < 0, then the 
tilted phase Lpi occurs, otherwise one has the non tilted 
L/3 phase below the main transition temperature. 

Keeping this in mind, we choose the value of \C\ for 
lipid A (DPPE) to be less than \C Lp \ and that for lipid B 
(DPPC) to be greater than \Cl p \- The parameter values 
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FIG. 5: Calculated mean-field phase diagram of DPPC- DPPE 
mixture as a function of composition and temperature T. cf> is 
the mole fraction of DPPE. C DPP c = 4.66xlO _19 J, C DPP e 
= 3.58X10" 21 J, T* DPPC = -17°C, T* DPPE = -3°C, and J = 
4.0 k p T* . All other parameter values are as in fig[2] 




FIG. 6: Calculated mean-field phase diagram of a binary lipid mix- 
ture as a function of composition and temperature at J = 5.0k p T* . 
All other parameter values are as in fig[5] 



are chosen so as to reproduce the transition temperatures 
of DPPC and DPPE. T4 is chosen to be negative for 
DPPC and positive for DPPE. 

Using the above values of the parameters the binary 
phase diagram for DPPC-DPPE mixtures were calcu- 
lated, and the results are shown in fig|5] and fig|6] In 
figj5] the interaction parameter J = 4.0 fc^T* while in 
figE]J = 5.0fc s T*. 

The calculated phase diagram of figE] is found to be 
in good quantitative agreement with the experimental 
one for DPPC/DPPE mixtures [48| . and is very sim- 
ilar to what is obtained for other lipid mixtures such 
as, DMPC/ dipalmitoyl phosphatidylserine (DPPS) and 
DPPC/DPPS [8]. At higher temperatures, we have the 
L a phase for the entire range of composition. With de- 
creasing temperature the Ppi phase appears very close 



to the pure DPPC axis, and is bounded by coexistence 
regions with the Lpi, Lp and the L a phases. The Ppi - 
L a coexistence region is very narrow on the scale of the 
figure and collapses almost into a line. At low tempera- 
tures there is a continuous Lpi - Lp transition at <f> ~ 0.2, 
which reflects the change in the sign of T4 according to 
eqnfT4l 

Fig. |5] shows the phase diagram for a binary lipid 
system in which the constituent lipids show a very 
high tendency to segregate and resembles the experi- 
mentally obtained one for dielaidoyl phosphatidylcholine 
(DEPC)/DPPE mixtures J4{|. The two lipids are com- 
pletely miscible only in the L a phase. The main qualita- 
tive difference from the previous case is the large region 
of Lpi - Lp coexistence at lower temperatures. 



2. Tilt-Tilt Mixtures 

In this section, we consider the case where both the 
constituents exhibit the ripple phase between the L a and 
Lpi phases. This would correspond to mixtures of lipids 
such as DPPC and DMPC. We, therefore, choose the 
value of C for the two lipids to be larger than Cl p and the 
sign of T4 to be negative for both of them. The reference 
temperature for A and B is chosen to be, T\ = —35° C 
and Tg = —17° C, respectively. The value of T4 is taken 
to be the same for both the lipids. Using these values of 
the parameters the binary phase diagram was calculated, 
and the results are shown in fig 7(a) fig |7(b)| and fig 7(c) 
for different values of J. 

When J = 3.0k B T* (fig 7(a) ) the three phases extend 
throughout the entire range of composition. These re- 



gions are separated by two coexistence regions of L a - 
Ppt and Pfji - Lp/ phases. As J is increases to 4.0 UbT* 
(fig |7(b)[ ), these coexistence regions almost touch each 
other, confining the P$t phase to a narrow strip. With 
further increase in the value of J (= 5.0 ksT*) two-phase 
coexistence regions dominate the phase diagram with the 
pure phases confined to compositions close to the two 
pure lipids. For this value of J the two lipids show a 
miscibility gap in the L a phase, with the appearance of 
a critical point. As discussed in the next section, the 
calculated phase behavior at low values of J closely re- 
sembles that observed in mixtures of very similar lipids. 
Although the complex phase diagrams calculated for high 
values of J have not yet been reported in experiments, 
we believe that such phase behavior should be seen in 
binary mixtures of very dissimilar lipids, such as those 
with hydrocarbon and fluorocarbon chains. 



IV. PHASE BEHAVIOR OF DPPC-DMPC 
MIXTURES 

Small angle scattering techniques can in principle de- 
tect microscopic phase separation in the plane of bilayers, 
if there is sufficient contrast in the scattering densities of 
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the two phases. However, even in the absence of such 
contrast, macroscopic phase separation can easily be de- 
tected from non overlapping reflections in the diffraction 
pattern coming from the individual phases. On the basis 
of the diffraction patterns we have determined the partial 
phase diagram of the binary mixture consisting of DPPC 
and DMPC. The use of oriented samples has an advan- 
tage that the in-plane ordering of molecules can be easily 
inferred from the wide angle reflections. For example, in 
the gel phase of both DPPC and DMPC, one observes 
two wide angle reflections, one on-axis (q z = 0) and the 
other off-axis {q z ^ 0), coming from the quasi hexagonal 
lattice of hydrocarbon chains of the lipid molecules. The 
positions of these peaks indicate that the molecules are 
tilted with respect to the bilayer normal and that the 
direction of the tilt is toward nearest neighbor [50| . 
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FIG. 7: Calculated mean-field phase diagram of a binary lipid 
mixture as a function of their relative composition </> and Temper- 
ature T. Both the lipids constituting the mixture exhibit a low 
temperature tilted phase. The lipid interaction parameter is set to 
bc[(a7lJ = 3.0k B T*,\(b)\j = 4.0k B T*,\(d\j = 5.0k B T*. C DMPC 
— f; onvin-19T QM ,i r>~„„~ — a c\c\^ in— i9 t All ^fV.^». r™».o mQ + Q ,. 



5.20xlCr ia J and C DPP c 
values are as in fig [2] 



4.66x1(5 J. All other parameter 



A. Materials and Methods 

l,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) 
and 1 ,2-dimyristoyl-sn-glycero-3-phosphocholine 

(DMPC) were purchased from Avanti Polar, and 
were used as received. Aligned multibilayers of binary 
lipid mixtures were prepared as follows: A concentrated 
solution of DPPC-DMPC binary mixtures (dissolved 
in a mixture of chloroform and methanol (1:1)) in 
designated molar ratios was deposited on the outer 
surface of a clean cylindrical glass substrate (radius of 
curvature of 9 mm). After deposition, the samples were 
placed overnight under vacuum to remove any remaining 
traces of the solvent. Subsequently, they were kept in 
a water-saturated atmosphere and were hydrated for 
a couple of days to obtain a stack of bilayers oriented 
parallel to the surface. 

Cu K a (A = 1.54 A) radiation from a rotating an- 
ode x-ray generator (Rigaku, Ultra X18) operating at 48 
kV and 70 mA and rendered monochromatic by a mul- 
tilayer mirror (Xenocs) was used to illuminate the hy- 
drated sample kept inside a sealed chamber with two my- 
lar windows. The chamber temperature was controlled 
using a circulating water bath to an accuracy of ± 0.1 °C 
and the relative humidity (RH) inside it was maintained 
at 98 ± 2%, by keeping a reservoir of water. The axis of 
the cylindrical substrate was oriented perpendicular to 
the incoming x-ray beam, such that the beam is incident 
tangentially to the sample. 

Diffraction patterns were recorded on a 2D image plate 
detector of 235 mm diameter and 0.1 mm pixel size (Mar- 
research) . All the samples were first heated to a temper- 
ature above the main transition temperature of DPPC 
and the diffraction patterns were recorded during cool- 
ing from the L Q phase. The sample temperature and the 
RH close to the sample were measured with a thermo- 
hygrometer (Testo 610) inserted into the chamber. 




FIG. 8: Small-angle diffraction patterns of aligned samples of 
DPPC - DMPC mixtures in the L a phase at T = 45.0°C, RH = 
99.9%, <t> = 0.5 (left), P , phase at T = 30.1°C, RH = 99.9%, 
(j> = 0.5 (right). The bilayer normal is vertical. 




FIG. 9: Diffraction pattern of an aligned sample of DPPC - DMPC 
mixtures in the Lpi phases, showing both the small-angle lamellar 
peaks due to the periodic stacking of the bilayers and the wide- 
angle peaks due to the in-plane ordering of the tilted chains in the 
bilayer. T = 14.9°C and RH = 99.9%, (j> = 0.5 . The bilayer 
normal is vertical. 



B. Results and discussion 



lattice. The wavelength of the bilayer height modula- 
tions in this phase is typically ~ 150 nm. The lower 
temperature Lpi phase shows more lamellar reflections 
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FIG. 10: Phase diagram of DPPC-DMPC bilayers at 98% RH 
determined from x-ray diffraction data, y, • and A correspond to 
the L a , Pa/ and Lai phases, respectively. 



in the small angle region compared to the L a phase due 
to decreased thermal undulations of the rigid bilayers. In 
addition sharp wide angle reflections are observed in this 
phase due to the ordering of the chains in the plane of 
the bilayer. The presence of a on-axis [q z = 0) and an 
off-axis (q z y^ 0) reflection show that the chains are tilted 
, with the tilt direction along the nearest neighbour. The 
DPPC-DMPC phase diagram closely resembles that cal- 
culated from our model for low values of J. The low value 
of J in this case is understandable since the two lipids 
are very similar with the only difference being that each 
chain of DPPC is longer by two CH3 groups. 



Figure [TO] shows the partial phase diagram of DPPC - 
DMPC mixtures at 98 ± 2% Rh. The three pure phases 
extend over the complete composition range, with two 
two-phase regions separating them. The phases were 
identified from their diffraction patterns. The L a phase 
gives a few reflections in the small angle regions, which 
correspond to different diffraction orders from a periodic 
lamellar stack. The lamellar periodicity in this phase is ~ 
5.6 nm for DPPC and ~ 5.0 nm for DMPC; the mixtures 
have intermediate values depending on the composition. 
The wide angle reflection due to chain ordering is very 
diffuse in this phase and cannot be seen in the present 
geometry for exposure times of the order of 10 min. Ad- 
ditional satellite reflections appear in the Ppi phase due 
to the periodic height modulations of the bilayer. These 
reflections can be indexed on a two-dimensional oblique 



V. CONCLUSION 

We have developed a continuum theory of phase transi- 
tions in single- and two-component achiral lipid bilayers. 
This model is found to reproduce all the salient struc- 
tural features of the ripple phase occurring between the 
fluid L a phase and the gel Lpi phase. In addition, it 
predicts the existence of a novel variant of this phase, 
which has not been experimentally observed so far. This 
theory has been used to calculate some generic phase di- 
agrams of binary phospholipid mixtures. We have also 
determined the phase diagram of DPPC-DMPC mixtures 
from x-ray diffraction studies on aligned multilayers. It 
is found to be in good agreement with that calculated 
from the model for weakly segregating species. 
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